Universal and nonuniversal contributions to block-block entanglement in 

many-fermion systems 
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We calculate the entanglement entropy of blocks of size x embedded in a larger system of size 
L, by means of a combination of analytical and numerical techniques. The complete entanglement 
entropy in this case is a sum of three terms. One is a universal x and L-dependent term, first 
predicted by Calabrese and Cardy, the second is a nonuniversal term arising from the thermody- 
namic limit, and the third is a finite size correction. We give an explicit expression for the second, 
nonuniversal, term for the one-dimensional Hubbard model, and numerically assess the importance 
of all three contributions by comparing to the entropy obtained from fully numerical diagonalization 
of the many-body Hamiltonian. We find that finite-size corrections are very small. The universal 
Calabrese-Cardy term is equally small for small blocks, but becomes larger for x > 1. In all inves- 
tigated situations, however, the by far dominating contribution is the nonuniversal term steming 
from the thermodynamic limit. 

PACS numbers: 03.67.Mn, 71.10.Fd, 03.65.Ud, 71.10.Pm 



I. INTRODUCTION 

Entanglement is one of the most studied and least in- 
tuitive features of quantum mechanics. Many aspects of 
it are still not fully understood. Part of the difficulty is 
that in itself entanglement is not an observable quantity. 
Rather, entanglement is a property of the quantum me- 
chanical state, defined with respect to some set of degrees 
of freedom. For different degrees of freedom, and differ- 
ent states, entanglement is characterized and quantified 
in different ways. 

For mixed states, described by a density operator, 
many alternative measures of entanglement are still un- 
der study. For pure states in a bipartite system, de- 
scribed by a wave function, on the other hand, a near- 
consensus has emerged that the entanglement entropy 
is a suitable entanglement measure. Having identified a 
suitable measure, the task at hand then changes to eval- 
uating it and analysing its behaviour in various physical 
systems, in order to extract information that can be use- 
ful in quantum information processsing and computing. 

The present paper is concerned with this task in the 
particular case of strongly interacting electrons in a 
finite-size chain. Our interest is in separating univer- 
sal and system-specific contributions to the entanglement 
entropy, quantifying their relative importance, and inves- 
tigating their behaviour as a function of system parame- 
ters. Specifically, we consider a quantum chain of length 
L divided in a subsystem A of size x, and a subsystem B 
of size L — x, and calculate the entanglement entropy [l| 

S(x,L) = -Tr[ PA \og 2 (p A )}, (1) 

where the reduced density matrix p A = Trs[p\ is ob- 



tained from the density matrix of the full system, p, by 
tracing over the degrees of freedom of subsystem B. For 
interacting many-particle systems the full density ma- 
trix is almost impossibly difficult to obtain. For suitable 
model Hamiltonians, and not too many particles, how- 
ever, fully numerical diagonalization is within reach, and 
can be used to calculate p A and S. Before embarking on 
such a numerical calculation for a specific system, how- 
ever, it is useful to recall general properties of S(x, L) 
that were uncovered in ground-breaking analytical work 
of Calabrese and Cardy [2| ■ 

These authors find that the entanglement entropy of a 
subsystem of size x, embedded in a larger gapless system 
of size L ^> x, consists of two distinct terms: a universal 
term depending only on x and L, and a nonuniversal 
term that depends on system-specific parameters, but is 
independent of x and L 0, H|- Analytical expressions 
for the universal term were obtained by Calabrese and 
Cardy (CC) 0, found to be in agreement with partial 
results obtained earlier in Refs. J3|,y,l5j, and were further 
analysed in, e.g., Refs. [1, 0, H, Ef- 

For periodic boundary conditions and L 3> x 3> 1, 
these authors find 



S(x,L) = -log 2 
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where c is the central charge (conformal anomaly) of the 
system and si is a nonuniversal term whose magnitude 
and dependence on system parameters remain undeter- 
mined in the approach of Refs. 0,11, 0,@. If the condition 
L 3> x 3> 1 is not satisfied there may be additional finite- 
size corrections, not contained in the CC analysis. While 
the identification of universal terms is one of the principle 
goals of statistical physics, any quantitative application 
to realistic models or to actual materials and devices, 
depends crucially on information about the nonuniversal 
terms. With a view on future realizations of quantum 
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computing and quantum information processing in sys- 
tems of interacting particles, we therefore now embark 
on the task to extract information about the nonuniver- 
sal function s\ and on possible finite-size corrections for 
realistic models of such systems. 

In this paper, we focus on the one-dimensional 
fcrmionic Hubbard model. For this model, we (i) nu- 
merically assess the magnitude of the universal Cardy- 
Calabrese (CC) term and the nonuniversal s% term for 
realistic values of system parameters, arriving at the un- 
expected conclusion that the universal term is only a 
small correction to the much larger nonuniversal term; 
(ii) obtain an analytical expression for s\ of the Hubbard 
model, allowing us to study its dependence on various 
system parameters; and (iii) compare the analytical ex- 
pression with numerical results obtained by full diago- 
nalization of the many-body Hamiltonian, finding satis- 
factory agreement, both for single-site and block-block 
entanglement, in interacting and noninteracting systems 
of various sizes and densities. 



II. UNIVERSAL CONTRIBUTION TO THE 
ENTROPY VERSUS EXACT ENTROPY 

In this section we numerically calculate the exact 
entanglement entropy of the one-dimensional finite-size 
Hubbard model and compare it to the analytical predic- 
tion made by keeping only the universal term in the CC 
formula. 



TABLE I: Universal contribution to the single-site (x — 1) 
entanglement of a chain of size L with periodic boundary 
conditions, compared to data from numerically exact full di- 
agonalization, at U = 4. The last column is the deviation, in 
percent, of the universal contribution from the exact numeri- 
cal value. 



n 


L 


guniv 


gcxact 


deviation (%) 




4 


-0.101 


1.541 


106.5 


0.5 


8 


-0.025 


1.564 


101.6 




12 


-0.011 


1.576 


100.7 




4 


-0.050 


1.594 


103.1 


1.0 


8 


-0.012 


1.701 


100.7 




12 


-0.005 


1.718 


100.3 



TABLE II: Universal contribution to the block-block entan- 
glement of a chain of size L = 10 with U = and n = 1, 
compared to values for the full entropy extracted for the same 
system from Ref. [l5|]. The last column is the deviation, in 
percent, of the universal contribution from the reference value. 



5 m 



,5" : ' 



deviation (%) 



-0.016 
0.602 
0.910 
1.065 
1.114 



2.00 
2.69 
3.02 
3.17 
3.22 



100.8 
77.6 
69.9 
66.4 
65.4 



S univ {x,L) 



log 2 
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(3) 



The one-dimensional Hubbard model, one of the most 
widely used models of strongly interacting particles [ToL 
[TTI . Il2| , is described by the Hamiltonian 



H 



(4) 



where t is the hopping between neighbouring sites, U 
is the on-site particle-particle interaction and cJ CT and 
Ci a are (fermionic) creation and annihilation operators 
of particles at site i with spin a. The system described 
by this Hamiltonian is completely characterized by its 
size L, interaction U [H| and either the particle num- 
ber N < 2L or the particle density n = N/L [l4j. For 
small L, this Hamiltonian can be diagonalized numeri- 
cally. Since this involves no approximation other than 
the use of finite-precision numbers on a computer, we 
follow common terminology and denote this as exact di- 
agonalization. The resulting eigenfunctions can be used 
to construct the density matrix, and from this the entan- 
glement entropy can be extracted. 

A quantitative comparison between the resulting 5' exact 
and 5 unlv , as given by Eq. ([3]), is presented in Table HI for 
various different choices of system parameters and x = 1. 



Interestingly, the universal term makes only a very small, 
negative, contribution to the exact single-site entropy. 

We conclude from this analysis that for the single-site 
entanglement the universal CC term is hardly relevant 
quantitatively: when L is large enough for the full CC 
formula to become asymptotically exact (recall that it 
was derived for L x), the universal term is already a 
vanishingly small correction to the nonuniversal term. 

From Eq. ([3]) it is clear that the universal term in- 
creases as a function of block size x. The interesting 
question is then by how much it grows relative to the 
nonuniversal, system-specific, term s\. In order to in- 
vestigate this quantitatively, we fix L at 10 sites, and 
calculate S(x, L = 10) as a function of the size x of the 
subsystem B. Benchmark data for comparison at x > 1 
are extracted from Ref. [l5| , which deals with block-block 
entanglement in the extended Hubbard model. By set- 
ting the parameter V of that model equal to zero, ref- 
erence values for S(x, L = 10) can be extracted from 
Figure 3 of that work. Below we do refer to these data as 
"exact" because they were also obtained by numerically 
exact full diagonalization, but we note that we extracted 
them graphically from Figure 3 of Ref. [HI]. The differ- 
ence between these "exact" values and S nnlv , which is the 
effect we are after, is clearly much larger than any pos- 
sible error of the benchmark data, which can therefore 
safely be used for comparison. 



3 



as 

ex o.o 





1 










- 




s 

o 
















E 














fi o 












H 


o 












s 















X 



FIG. 1: Derivative of the block-block entanglement entropy 
with respect to block size, obtained analytically from Eq. ([2]) 
[open circles] , numerically from five data points obtained from 
Eq. ((2} [open squares], and numerically from the five data 
points in Table ITT1 



diction of the trend as a function of x and the numerical 
results, both obtained from five data points. Evidently, 
the behaviour of the CC expression and of the numeri- 
cal data is very similar. The small differences remaining 
between crosses and squares are due to finite-size correc- 
tions to the CC expression (more on these below) and the 
intrinsic error bar of the reference data from Table [TXJ 

Since the behaviour of the derivative dS/dx is very 
closely reproduced by the CC expression, we conclude 
that the much larger differences observed in Table [H] for 
the entropy itself, must be almost entirely due to the 
term S\. This analysis thus points to the importance of 
nonuniversal terms, which remain undetermined in the 
CC approach. 

Additionally, it should be noted that the CC formula 
was derived for large L ^> x, while in order to be able 
to compare to data from full numerical diagonalization 
we evaluate it for L < 13. Finite-size corrections to the 
CC formula are another possible explanation for the large 
fraction of the exact entropy not recovered by the uni- 
versal term only. In the next sections we attempt to 
disentangle and quantify these two distinct effects, by 
deriving an analytical expression for si of the Hubbard 
model, and comparing it to the same set of exact data. 



The data in Table |TT] show that for larger blocks the 
universal term makes a more noticeable contribution to 
the block-block entanglement entropy. While for x = 1 
the deviation of the universal term from the numerical 
value is very similar to that of TablelH it becomes smaller 
for x > 1. However, quantitatively, it is still amounts to 
less than half of the exact entropy. As an example, even 
for maximal subsystem half the size of the 

complete system, it is still only about one third of S exact . 
Part of this difference must be due to si. 

The si term can be eliminated by studying not 5, but 
its derivative as a function of block size x. This behaviour 
can be evaluated in various ways. First, we calculate an- 
alytically the derivative dS(x, L)/dx of the CC expres- 
sion ([2]), by treating x as a continuous variable. Results 
for integer values of x are represented by open circles in 
Fig. [1] Note that since si is taken to be a constant, this 
derivative samples only the universal term. Second, we 
calculate the numerical derivative of the five data points 
for x = 1,..,5 collected in Table iHl Results are repre- 
sented by crosses in Fig. [TJ A direct comparison between 
both sets of data is marred by the intrinsic inaccuracy 
of a numerical derivative. To minimize this problem, we 
also obtained the derivative of the CC formula © by 
evaluating that expression numerically at x — 1, ..,5 and 
taking the numerical derivative of the resulting values. 
This set of data, represented by open squares, is directly 
comparable to the numerical derivative of the benchmark 
data in Table HU 

The deviation between the open circles and open 
squares is thus between analytical and numerical 5-point 
derivatives of the same function. The difference between 
the open squares and the crosses is between the CC pre- 



III. ANALYTICAL EXPRESSION FOR THE 
BLOCK-ENTANGLEMENT ENTROPY 

In the thermodynamic limit, L — > oo, the CC expres- 
sion ([2]) reduces to 

S(x,L -v oo) = -log 2 (a;) + si, (5) 

from which the function si of the model under study 
can be determined once the entanglement entropy of this 
model is known in the thermodynamic limit. Since Si 
does not depend on x, we are free to evaluate it for any 
convenient value of x. In the special case of single-site 
entanglement (x — 1), the logarithm on the right-hand 
side vanishes, and 

si = S(x = 1, L — * oo), (6) 

which implies 

S(x,L -> oo) = ^ | c log 2 (a;) 
iS(l, L — » oo) 3 Si ' 

This identification neglects possible finite-size corrections 
arising from the fact that the CC expression was derived 
only for L»i> 1. The difference between numerical 
data and predictions of the preceding equation allows one 
to estimate the size of such corrections. 

Next, we apply this procedure to the Hubbard chain. 
Recent research has resulted in a complete physical pic- 
ture of and explicit expressions for S(x = 1,L — * oo) 
G3, EH ■ Specifically, the single-site entanglement en- 
tropy for the Hubbard model in the absence of external 
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electric or magnetic fields is given by [16|, [l7|, [IE 



it de \ 

S{x = 1, L -> oo; n, U) = -2 ( - - — J log 2 

de\ 

mi log2 



n de 
2~dU 
de 

n + W 



,(8) 



where the ground-state energy per site, e = Eq(u, U)/L, 
can be obtained from the Bethe-Ansatz integral equa- 
tions o, m d, d . 

By combining this Bethe-Ansatz based expression for 
S(x = 1, L — > oo) with the CC formula, based on confor- 
mal field theory, we obtain 



S(x, L; n,U) — — log. 



3 
de 
dU 



log 2 



1-n 



m ) log2 



L . / 7TX s 



n de 
2~dU 
de 
dU 



1 - n 



(9) 



which is the sum of a universal term, depending only on 
geometry {x and L), and a term depending on specific 
system parameters (n and U). This explicit expression 
allows us to investigate the actual size of each of these 
terms under realistic circumstances, by obtaining e(n, U) 
numerically from the Bethe-Ansatz integral equations, 
and evaluating Eq. ([9]) as a function of x, L, n and U. 



IV. FULL ENTROPY VERSUS EXACT 
ENTROPY 

Figure [2] contains a comparison of our numerical data 
for x — 1 with our analytical expression ([9|) for the spe- 
cific case of a Hubbard chain with U = 4 and n = 0.5. 
Exact data are given only for L = 4,8, 12, because odd 
particle numbers N would result in a finite magnetization 
(which is not included in Eq. ([5])) and L > 13 is already 
too large for full exact diagonalization on our comput- 
ing equipment. The available data, however, are clearly 
sufficient to conclude that trend and magnitude are the 
same for both analytical and numerical data. 

The quantitative deviation observed in Fig. [2] between 
analytical and numerical results for small L is due to the 
fact that the CC formula was derived for large L, whereas 
data from full numerical diagonalization are only avail- 
able for L < 13. For small L, there may be additional 
L-dependent terms in the full expression for the entropy, 
which go to zero as L increases. These are the finite- 
size effects, referred to above. The data in Fig. [5] show 
that already for L as small as 12, such possible small-L 
corrections are negligible. 

Figure [3] extends this analysis to larger block sizes, by 
evaluating Eq. ^ as a function of x, and comparing to 
the same set of reference data at x > 1 and U = 0, used 
in Sec. HIl The overall agreement between expression © 




FIG. 2: Single-site entanglement entropy as a function of sys- 
tem size L for the Hubbard model. Open squares: analytical 
results from our Eq. ©. Crosses: numerical data obtained 
by diagonalizing the many-body Hamiltonian. 



TABLE III: Same as Table U but with S univ (x,L) replaced 
by our expression ([9}, comprising S unlv (x,L) and expression 
© for si(n,U). 



n 


L 


S(x = 1,L) 


exact 


deviation (%) 




4 


1.475 


1.541 


4.3 


0.5 


8 


1.551 


1.564 


0.8 




12 


1.565 


1.576 


0.7 




4 


1.678 


1.594 


5.3 


1.0 


8 


1.716 


1.701 


0.9 




12 


1.723 


1.718 


0.3 



and the benchmark data is rather satisfactory. The re- 
maining deviations now have two distinct sources. One 
is, as above, the use of the CC expression at rather small 
L. The other is the imprecision in the extraction of the 
benchmark data from the figure presented in Ref. [l5l |. 
Still, there can hardly be any doubts that the depen- 
dence on block size x is reproduced correctly. Clearly, 
in Eq. © this dependence comes exclusively from the 
universal term, which makes a much more pronounced 
contribution for x > 1 than it made for x = 1. 

Tables lHll and llVl compare expression © to benchmark 
data at x = 1 and x > 1. The last column of Tables ITTT1 
and llVl shows that expression §§§ practically exhausts the 
exact entropy, both for single-site entanglement and for 
block-block entanglement. As before, we attribute the 
remaining small differences, of order ~ X%, to finite-size 
corrections, contained in the numerical data for small L 
but not in the CC expression derived for large L. Since 
for some values of n the full expression ^ predicts more 
than 100% of the exact entropy, these finite-size correc- 
tions must alternate their sign as a function of n. 
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FIG. 3: Block-block entanglement entropy as a function of the 
block size x for the noninteracting Hubbard model (L = 10, 
U = and n — 1). For x > L/2 we obtained S from the 
symmetry relation S(x,L) = S(L — x,L), and for x < L/2 is 
calculated from Eq. © ■ Benchmark data have been extracted 
from Ref. [l]|. 



TABLE IV: Same as Table HH but with S univ (x,L) replaced 
by our expression ©, comprising 5 unlv (a;,L) and expression 
© for si(n,U). 



X 


S(a;,L) 


exact 


deviation (%) 


1 


1.984 


2.00 


0.8 


2 


2.602 


2.69 


3.3 


3 


2.910 


3.02 


3.6 


4 


3.065 


3.17 


3.3 


5 


3.114 


3.22 


3.3 



As a second test, we have also fitted the exact data for 
L = 10 and x < 5 in Tables |TT] and IIVI with expression 
@, treating c and Si as fitting parameters. The result 
is c = 2.17 ± 0.02 and Si = 2.02 ± 0.007. Since in the 
situation of these Tables the exact values are known to 
be s\ = c = 2, this fit again illustrates the smallness of 
finite-size corrections to the CC formula <j2j) for x and L 
outside the range L>2!>1 ! where it becomes exact. 



V. CONCLUSIONS AND OUTLOOK 

We can summarize our combined analytical-numerical 
analysis in the statement that the block-block entangle- 
ment entropy for any system size L and block size x is 
given by 



-S(x 



S(x,L;n,U) = S unw (x,L) 
l,L-Kx>;n,U) + AS(x, L; n, U), (10) 



where the first term is the universal CC term of Eq. ([3]), 
depending on x and L only, and the second is the system- 
specific term s\, depending on n and U, which we ex- 
tract from the thermodynamic limit The third by 
definition comprises all possible finite-size corrections 
not resulting from the CC analysis and neither con- 
tained in the infinite-size limit leading to the identifi- 
cation si (n, U) = S(x = 1, L oo; n, U). 

By comparing this full expression to the data in Ta- 
bles U to IIVI and the trends visible in Figs. Q] to O we 
conclude that the nonuniversal term s\ = S(x = 1, L — * 
oo; n, U) quantitatively dominates the physics of the en- 
tanglement entropy of the Hubbard model in all inves- 
tigated situations. For x = 1, the universal term and 
the finite-size corrections are of comparable magnitude, 
0(1%), and essentially negligible relative to the nonuni- 
versal term, even for rather small L. For x > 1, the 
universal term is of order O(10%) of the full entropy, 
while finite-size corrections remain 0(1%). The trend 
as a function of x, on the other hand, is not at all af- 
fected by si, but dominated by the universal term, re- 
ceiving only small (~ 1 — 10%) finite-size corrections from 
AS(x,L;n,U). 

While all of this highlights the intellectual achievement 
of CC in having identified the universal contribution to 
the entanglement entropy, it also shows that if one wants 
to quantify the entropy in an actual material or device 
— a need that arises as soon as one considers using en- 
tanglement as a resource for quantum computing in real 
systems — a detailed description of system-specific fea- 
tures is unavoidable. 

One way to obtain an approximate system-specific de- 
scription also in systems with inequivalent sites is the 
local-density approximation (LDA) to density-functional 
theory (DFT), which locally applies results obtained in 
a spatially homogeneous system (with constant density 
n) in order to simulate the corresponding inhomogeneous 
system (with spatially varying density m ) . This approx- 
imation is commonly applied in ab initio calculations of 
the electronic structure of atoms, molecules and solids 
|20l fell |22| . |23j , where the local approximation (or one of 
its refinements) is made for the exchange-correlation en- 
ergy. In earlier work we proposed to apply the same strat- 
egy also to the Hubbard model [24[ , and suggested a sim- 
ple local-density approximation for the single-site entan- 
glement entropy [25|]. (For related applications of DFT 
concepts to the study of entanglement, see Refs. [26lfe"7j.) 

That approximation was constructed specifically for 
i=l. The present equation (flT))) . valid for all x, suggests 
two ways to extend the validity of the local-density ap- 
proximation to block-block entanglement. One is to sim- 
ply add to the LDA of Ref. H the term MS univ (x, L), 
which the original entropy LDA did not contain. Here 
M = L/x is the number of blocks in the system. The 
other is to take Eq. ifTO]) as entropy of the homogeneous 
reference system on which the local approximation is 
based. This leads to 



S 



LDA 



x,L;m,U] = MS UDiv (x,L 



6 



-S LDA [x = l,L;m,U], (11) 

x 

which differs from the ad hoc correction in the in- 
dependence of the first term and the x-dependence of 
the second. Future work directed at block-block entan- 
glement in spatially inhomogeneous systems is expected 
to identify which of these two extensions is more reliable. 
Interesting inhomogeneities, in this context, include im- 
purities, defects, spatial modulations of the system pa- 
rameters, confining potentials, etc. 



A combination of the methodologies employed in 
Ref. [25[ and in the present paper thus allows to analyse 
and quantify the entanglement entropy in a wide variety 
of spatially inhomogeneous many-body systems. Such 
analysis, and an extension of these investigations to spin- 
polarized systems, is subject of future work. 
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